Realistic modeling of the electronic structure and the effect of correlations for 

Sn/Si(lll) and Sn/Ge(lll) surfaces 
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The correlated electronic structure of the submonolayer surface systems Sn/Si(lll) and 
Sn/Ge(lll) is investigated by density-functional theory (DFT) and its combination with explicit 
many-body methods. Namely, the dynamical mean-field theory and the slave-boson mean-field the- 
ory are utilized for the study of the intriguing interplay between structure, bonding and electronic 
correlation. In this respect, explicit low-energy one- and four(sp 2 -like)-band models are derived 
using maximally-localized Wannier(-like) functions. In view of the possible low-dimensional mag- 
netism in the Sn submonolayers we compare different types of magnetic orders and indeed find a 
120° antiferromagnetic ordering to be stable in the ground state. With single-site methods and 
cellular-cluster extensions the influence of a finite Hubbard V on the surface states in a planar and 
a reconstructed structural geometry is furthermore elaborated. 

PACS numbers: 73.20.At, 71.15.Mb, 71.10.Fd, 71.30. +h 



I. INTRODUCTION 



Strong electronic correlations lead to many of the 
most interesting phenomena in modern condensed matter 
physics, such as superconductivity, the metal- insulator 
transition or local moment magnetism. It is well known 
that the effect of these many-body effects depend heav- 
ily on the interplay between structure, bonding, and the 
degree of order within the effective dimensionality of the 
given problem. While the study of realistic bulk sys- 
tems in this context has been rather extensive in recent 
years, those investigations may only provide insight into 
the effect of dimensionality via reasonably justified quasi- 
dimensions due to large crystalline anisotropics. In this 
respect the study of surface systems has the advantage 
that by an intelligible tuning of the coupling between 
the adsorbate and the substrate, quasi-lowdimensionality 
can be created more efficiently. Besides, there are many 
surface-sensitive experimental probe techniques in order 
to reveal the electronic structure more directly than in 
the bulk. Hence the research on specific surface materials 
problems may add substantially to the understanding of 
the general problem of strong quantum correlations. 

There has been an increased interest in adsorbate 
systems involving semiconducting substrates since more 
than a decade^—. These often exhibit dangling-bond- 
derived surface states with rather small bandwidths. 
While Mott criticality is hardly seen on the free sur- 
faces because of structural reconstructions, certain adsor- 
bate atoms forming (sub)monolayers stabilize hybridized 
narrow-band surface states down to low temperatures. 
The most prominent of this kind are the so-called a-phasc 
surfaces^ where the canonical structure is described by a 
y/3x^/3R30° triangular array of adsorbate atoms within 
an 1/3 monolayer coverage on a (111) semiconductor sub- 
strate, resulting in a half-filled surface band. Due to the 
rather large interatomic distances between the adsorbate 
atoms (~7A) the latter becomes indeed rather narrow. 
Mott-insulating phases in such systems are believed to 



be realized in e.g., K/Si(lll)2, whereas a transition to a 
correlation-driven charge-density-wave phase takes place 
in Pb/Ge(lll)i. 

In this respect, the a-phase systems Sn/Si(lll) and 
Sn/Ge(lll) are of central interest due to the unusual 
properties they exhibit at low temperatures. It has been 
recently verified by Modcsti et &\X that the planar a- 
Sn/Si(lll) surface exhibits a metal-insulator transition 
below 60 K. No structural reconstruction appears to ac- 
company this transition, since it was observed^ that the 
s/3x V3i?30° periodicity is stable at least down to around 
6 K. The question of magnetic ordering within the insu- 
lating regime, i.e., the formation of an antiferromagnetic 
(AFM) Mott-insulating state, has been raised because of 
band foldings revealed by low T photoemission experi- 
ments^^ leading to an 3x3 periodicity. In contrast, the 
structurally and electronically very similar a-Sn/Gc(lll) 
surface shows vastly different behaviour. There a transi- 
tion from the v / 3xv / 3-R30° phase at room temperature to 
an 3x3 symmetry below 200K takes place. Two compet- 
ing ground-state configurations of the Sn atoms seem to 
exist in the latter temperature ranged, with the so called 
2D-1U (i.e., two Sn atoms down, one Sn atom up with 
respect to the planar triangular structure) state in favor 
of 6 meV/adatom against the competing 1D-2U configu- 
ration^. The -\/3x v / 3i?30° periodicity observed at room 
temperature is understood as a result of the rapid fluctu- 
ations of the system between these two states^. It is still 
a matter of debate whether or not the system displays a 
surface Mott transition similar to the Sn/Si(lll) system. 
The finding of a Mott insulating phase in a-Sn/Ge(lll) 
below 20K by Cortes et al*£ has apparently not been 
confirmed by other group a 11 ' 12 . 

The electronic structure of the a surface phases 
poses an interesting problem in the context of strong 
electron correlation. Early phase-diagram studies of 
such adlaycr structures within the Hartree-Fock ap- 
proximation showed the possibility for various order- 
ing^. Since the narrow-band surface state is mainly 
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derived from a hybridization of the Sn(5p z ) state with 
the underlying Si states, one deals with Coulomb cor- 
relations in an effective 5p system (assuming a local 
Hubbard- like interaction). Calculations based on stan- 
dard density-functional theory (DFT) sho w 14 ' 15 in the 
case of ^/3x\/3R30° Sn/Si(lll) indeed an isolated half- 
filled surface band of width cV. Additionally, 
the named 3x3 reconstruction in Sn/Ge(lll) is verified 
within DFT— . Correlation effects beyond the local 
density approximation (LDA) and the generalized gra- 
dient approximation (GGA) have been investigated by 
Profcta and Tosatti2£ within the LDA+U method^ 1 -, es- 
tablishing a Hubbard U for the single Sn(5p z ) orbital. In 
this scheme, Mott-insulating Sn/Si(lll) (with assuming 
ferro-/ferrimagnctic order) is reached for [7~2 eV, but 
the authors find a £/~4 eV from constrained calculations 
more appropriate, also to cope with a gap size of order 
~0.3 eV. This rather large value for the local Coulomb in- 
teraction differs from [/ e ff— 1.15 eV obtained in elder con- 
strained LDA calculations designed for a minimal model 
describing the surface bancU^. Concerning a-Sn/Si(lll) 
there are speculations^ ' 13 ' 20 about a realization of the 
quasi-twodimcnsional (2D) correlated triangular lattice 
problem, including the possiblity of antiferromagnetism 
or spin-liquid physics and d-wave superconductivity. 

In this work we want to investigate the importance of 
electronic correlations and their interplay with the struc- 
tural data in the ct-Sn/ (Si,Ge)(lll) phases. By means of 
a combination of realistic DFT band-structure schemes 
with advanced many-body techniques, the aim is to clar- 
ify the differences stemming from the (Si,Ge)(lll) sub- 
strates and to reveal to which extent Coulomb correla- 
tions can give rise to the observed and perhaps still-to- 
be obeserved phenomena. Albeit some relevant work 
has already been performed in this direction, namely 
Ref. ll3lfT5ll20l . there remain still many open questions. 
For instance, the minimal Hubbard model derived by Flo- 
res et alr^ has not truly been numerically treated and the 
arguments given concerning the influence of structural re- 
constructions on a possible Mott-criticality lack a local 
picture of the involved orbitals. The LDA+U method 
utilized in Ref. [2^ is designed for long-range-ordered insu- 
lating states and neglects as a static technique quantum 
fluctuations (even in the Mott state) . Hence correlations 
in the metallic regime are usually described incorrectly 
and magnetic tendencies are often overestimated. More- 
over, the experimental data is far from being conclusive, 
e.g., concerning the appearance of local moment physics 
and eventual magnetic ordering. 



II. THEORETICAL APPROACH 

Realistic theoretical schemes for correlated condensed 
matter, combining traditional band-structure approaches 
with explicit many-body techniques have been quite suc- 
cessful in the last decade in dealing with various prob- 
lems in strongly correlated physics. The most promi- 



nent of such schemes is the combination of DFT with 
the dynamical mean-field theory (DMFT), the so-called 
LDA+DMFT 2 ^ 3 . framework. While in standard DFT 
correlation effects are only taken into account in an av- 
eraged way by making reference to a homogeneous (-like) 
electron gas, the named combined approach allows for 
explicit many-body effects on an operator level by still 
keeping important band-structure details from LDA. The 
DMFT techniqu o 24 ' 25 is able to incorporate all onsite 
quantum fluctuations and thereby may describe quasi- 
particle (QP) as well as atomic excitations (i.e., Hub- 
bard bands) on an equal footing. A alternative com- 
bined approach is given by interfacing the Gutzwiller- or 
the slave-boson mean-field technique with DFT— ~— . In 
this more simplified treatment the QP lifetime remains 
infinite, thus omitting the full frequency dependence of 
the self-energy. Hence only low-energy features may be 
addressed in the spectral function and hence high-energy 
Hubbard bands are not accessible. However importantly, 
the local atomic multiplets are still present with an ef- 
fective static character in the generalized theory 2 ^— . 

Due to the subtle dependencies of the oc- 
Sn/(Si,Ge)(lll) surface electronic structure on the 
structural facts, care must be taken in its proper 
determination. We used an implementation 3 ^ of the 
highly-accurate mixed-basis pseudopotential (MBPP) 
technique^ 3 - for this task. This DFT band-structure 
code employs normconserving pseudopotentials 3 ^ and 
an efficient combined basis consisting of plane waves and 
a few localized orbitals. Since the magnetism of these 
surface systems is very subtle, we additonally performed 
computations within the projector-augmented-wave 
(PAW) method 3 ^ for the specific study of magnetic 
ordering. Thereby one is able to lift the possible limi- 
tations due to the use of pseudo-crystal wave functions. 
Two implementations of the PAW formalism, namely the 
CP-PAW— code and the Vienna Ab-initio Simulation 
Package (VASP) 3 ^ code, were applied, which also allow 
for the investigation of noncollinear spin orderings. 
The actual calculations were performed by employing 
a slab geometry where care was taken in using a well 
converged lateral k-point mesh (up to a 25x25 grid for 
the magnetic structures). 

We combined the DFT approach with two many-body 
techniques, namely the DMFT with a Hirsch-Fye quan- 
tum Monte-Carlo impurity solver— and the rotationally 
invariant slave boson (RISB) metho d 30 ' 38 in saddle-point 
approximation. Note that in the single-site problem the 
RISB technique may also be understood as a quasiparti- 
cle impurity solution to DMFT. The employed mean-field 
version of the RISB method is in its character and na- 
ture of approximation very similar to the state-of-the-art 
Gutzwiller techniqu o 29 ' 31 ' 39 . 

For the actual interfacing of DFT with the many-body 
approaches a suitable downfolding procedure of the rele- 
vant problem to a local basis (the correlated subspace) is 
needed (see e.g., Ref. l40l - |42h . The choice is in most cases 
a matter of convenience, since many physically sound 
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frameworks that provide a representation of the band- 
structure within a minimal local Wannier(-like) basis are 
applicable. The coherent connection of the crystal prob- 
lem to such a tailored basis is the key point, rather than 
the pecularitics of the (restricted) projected local view- 
point itself. It is important to realize that the extended 
(crystal, surface or chain) problem, although composed 
of atoms, is not an atomic problem. Of course, in a sub- 
sequent many-body treatment the minimal interacting 
hamiltonian has to be adjusted to the characteristics of 
the chosen local basis. Up to now, the results of such 
many-body approaches are not exceedingly sensitive to 
the very details of the selected correlated subspacc, once 
the number and character of orbitals is agreed on. In 
that sense, the here utilized maximally-localized Wan- 
nier (MLWF) schem o 43 i 44 thus provides a reliable local 
basis for the materials under consideration. 



III. DFT INVESTIGATION AND MLWF 
DESCRIPTION 

A. Single-site Sn unit cell 

In order to investigate the electronic structure of 
the 1/3 monolayer of tin atoms on the semiconductor 
y/3x\/3R30° surfaces we first utilized slab geometries in- 
corporating a single Sn atom within the full unit cell. 
Thus the surfaces are modeled by a superccll with three 
bilayers of (Si,Ge) separated by sufficiently large vacuum 
regions. Note that in this section we model the Sn sub- 
monolayer in each case as flat, although a distortion in 
z direction is observed for true Sn/Ge(lll) (see section 
IIIIBI) . The Sn atoms are placed in the T\ sites of the 
surface while the bottom of the slab is saturated with hy- 
drogen atoms (see Fig.[TJ. Structural relaxations by min- 
imizing the atomic forces have been performed with fixed 




FIG. 1: (Color online) Left: Top view onto the employed 
V3x\/3R30° (red) and 3x3 (green) surface unit cells . Right: 
Side view of the slab geometry for Sn/Si(lll) surface (bot- 
tom) and the distorted 1U-2D Sn/Ge(lll) surface (top). 



position of the lowest substrate layer. The bulk chosen 
lattice constants and relaxed Sn-Sn nearest-neighbor dis- 
tances are (5.43, 6.65)A and (5.65, 6.93)A for Sn/Si(lll) 
and Sn/Ge(lll), respectively. An energy cutoff of 16 
Ryd for the plane-wave part of the mixed-basis set was 
used for all calculations. Localized functions were intro- 
duced for Si(3s,3p), Ge(3<i,4s,4p) as well as Sn(4d,5s,5p). 
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FIG. 2: (Color online) GGA-PBE DOS for Sn/(Si,Ge)(lll) 
m the V3xV3R30° structure, (a) Total DOS together with 
local Sn DOS and (b) shown via contributions from different 
atoms within the supercell (for the labeling of the substrate 
atoms see Fig.JS]). (c) Angular- momentum resolved plot of the 
Sn DOS. Note that in the latter plot the 5p x and 5p y curves 
lay on top of each other. 
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FIG. 3: (Color online) Top: Surface band structure for the 
Sn/Si(lll)-\/3x^i?30 o system, with fatbands for the Sn(5p) 
orbitals (blue/grey: p x ,p y ; red/lightgrey: p z ). The bulk Si 
band structure is shown in the background. Middle: De- 
rived Wannier band (cyan/lightgrey) for the 1-band model. 
Bottom: Derived Wannier-like bands (cyan/lightgrey) for the 
four-band model. 

Note that hence the semicore d electrons of Ge and Sn are 
treated as valence in our computations and relativistic ef- 
fects are included via the scalar-relativistic normconscrv- 
ing pscudpotentials. Since usually more appropriate for 
surface studies than LDA, we employed the GGA in the 
form of the PBE functional 4 ^ to the exchange-correlation 
term in the density functional. 

The density of states (DOS) shown in Fig. [2] displays 
for both surface systems a prominent structure close to 
the Fermi level of width 0.4 eV (slightly larger for 
Sn/Ge(lll)). Two subpeaks arc visible, one nearby the 
Fermi energy and the other at the upper energy edge in 
the unoccupied region. The major contributions to this 
low energy part stems from the orbitals of the Sn adatom 
hybridizing with the orbitals of neighbouring substrate 
atoms of the first and second layer (see Fig. ODd) . It be- 
comes further obvious from the angular-momentum re- 
solved DOS in Fig. that concerning the tin part the 
Sn(5p 2 ) orbital is dominantly responsible for the low- 
energy weight. Minor contribution to this is also added 
by the Sn(5s) orbital. 

As seen from the band-structure plots in Fig. [3] and 
Fig. 21 the low-energy DOS is associated with a single 
surface band within the gap of the bulk system. The so- 
called fatband resolution for the Sn(5p) orbitals, i.e., the 
identification of the contribution of a given orbital char- 
acter to a Kohn-Sham band at each k point via the width 
of a broadened coating 'band', underlines the dominating 
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FIG. 4: (Color online) As Fig. 02 but here for the (flat) 
Sn/Ge(lll)-V3x V3R30° system. 

p z weight on this half-filled surface band. A straightfor- 
ward MLWF construction may be applied to the single 
surface bands for both systems. The real-space Wan- 
nier function (see Fig. [5]) shows not only the p z -like lobe 
but displays additionally the threefold bonding aspects 
to the substrate atoms. A spread of 14.2 A 2 for this 
MLWF was obtained in the case of Sn/Si(lll). Albeit a 
rather similar low-energy state may be identified in the 
flat Sn/Ge(lll) system, it has to be noted that conven- 
tional DFT has its problems in describing germanium 
systems. For the bulk material, there is a nearly vanish- 
ing band gap within DFT— and additionally relativistic 
effects on the band structure are also importan t 46 ' 47 . Still 
in our PBE-GGA study a single Wannier band may also 
be extracted in this case (see Fig. [4} , however contrary 
to Sn/Si(lll) that band touches occupied levels at the T 
point. 

Due to the substrate geometry and the apparent Sn(5s) 
contribution at the Fermi level, an extended low-energy 
modeling based on sp 2 hybridized orbitals seems even 
more adequate for these systems. Indeed as shown in 
Figs. (CUE]), the effective bands from the corresponding 
four-band MLWF construction fit well to the full Sn(5p) 
fatband dispersion (Sn(5s) has nearly exclusively weight 
on the band at the Fermi level). The chemically more 
appealing sp 2 +p z viewpoint yields three of the four or- 
bitals having a dominant in-plane orientation with an 
120° angle between them, while a fourth one is now di- 
rectly reminiscent of the p z orbital (see Fig. [6]). This 
latter Wannier-like orbital has indeed major weight on 
the low-energy band close to £f- Compared to the one- 
band case, the spread of the Wannier-like functions are 
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FIG. 5: (Color online) Wannier orbital corresponding to the 
one-band model. The numbers are references for the PBE- 
GGA DOS plot in Fig. [2] 




FIG. 6: (Color online) Left: Wannier-like orbitals of the four- 
band model, (a-c) sp 2 -like orbitals and (d) p z -like orbital. 
Right: Top view on the spatial orientations of the three sp - 
like orbitals relative to the y/3x y/3R30" unit cell. 



now given as 13.6 A 2 for the sp 2 -like orbitals and 20.8 
A 2 for the remaining p z -like orbital. 

The hopping intcrgrals for both Wannier construc- 
tions are provided in Tab. |TJ In the one-band model 
the nearest-neighbor hopping with an absolute value 
|t|~45 meV is negative in accordance with the hole-like 
dispersion from Fig. ([3J [JJ around the T point. For 
the Sn/Ge(lll) system the more distant hoppings are 
slightly larger, whereas \t\ is somewhat smaller compared 
to the Sn/Si(lll) case. The onsite energies in the four- 
band model yield a crystal-held splitting between the sp 2 - 
like orbitals and the p z -like one of the order of 1.6 (0.9) 
eV for Sn/(Si,Ge)(lll). Though some caution has to be 
taken because of the named DFT-Ge problems, this illus- 
trates the more isolated p z -like level for the Si substrate 
and underlines the stronger hybridzation in the case of 
the germanium substrate. 

While the minimal one-band model is easily justified, 
the eligibility of the four-band model needs some ad- 
ditional comments. Albeit from the bonding character 




FIG. 7: (Color online) Illustration of the Wannier-like orbitals 
in the four-band model for the assignement of the hoppings 
in Tab. [I] The sketched orbitals 1, 2 and 3 correspond to the 
top, middle and bottom orbitals in the right part of Fig. [5] 
Orbital number 4 denotes the remaining p z -like level. 



direction 


100 


110 


200 


210 


Sn/Si(lll) 


44.6 


-18.4 


6.7 


-0.8 


Sn/Ge(lll) 


43.2 


-23.7 


7.3 


-1.8 



direction 


Sn/Si hoppings [meV] 


Sn/Go hoppings [meV] 
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-433.1 
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-358.3 
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-433.1 


546.4 


-358.3 




452.0 


452.0 


-634.8 


-318.6 


546.4 


546.4 


-433.1 


-358.3 




-318.6 


-318.6 


-318.6 


990.3 


-358.3 


-358.3 
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463.2 




-42.6 


-29.9 


29.4 


-26.7 


-61.5 


-37.8 


3.6 


-20.6 


100 


-29.9 


-42.6 


29.4 


-26.7 


-37.8 


-61.5 


3.0 


-20.6 




74.2 


74.2 


50.8 


-33.1 


33.8 


33.8 


2.9 


70.9 




85.0 


85.0 


-49.1 


16.7 


54.9 


54.9 


1.6 


29.2 




-0.7 


-5.9 


-8.4 


1.5 


-15.2 


-19.0 


-22.8 


1.3 


110 


12.1 


-9.5 


-5.9 


11.3 


2.9 


-22.8 


-19.0 


14.3 




3.8 


12.1 


-0.7 


-24.5 


-10.9 


2.9 


-15.2 


-4.6 




-24.5 


11.3 


1.5 


-9.4 


-4.6 


14.3 


1.3 


-10.3 




-1.2 


-1.4 


4.2 


-1.2 


3.6 


-4.2 


2.8 


3.6 


200 


-1.4 


-1.2 


4.2 


-1.2 


-4.2 


3.0 


2.8 


3.6 




-1.6 


-1.6 


8.0 


1.9 


-2.5 


-2.5 


0.6 


7.3 




5.3 


5.3 


-2.8 


-1.1 


12.6 


12.6 


6.4 


-10.5 




0.7 


0.3 


0.9 


-3.3 


-4.3 


-0.8 


-2.5 


0.8 


210 


2.2 


0.9 


1.3 


-0.5 


-1.9 


0.2 


-2.7 


-1.1 




-0.4 


-2.9 


0.0 


-2.6 


-2.4 


-3.7 


-0.9 


-0.7 
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-3.4 


1.4 


3.3 


1.8 


-3.0 


1.5 


2.9 



TABLE I: Hopping integrals up to the fourth-nearest neigh- 
bours for the minimal one-band case (top) and the 4-band 
(sp 2 +p z ) model (bottom). For the latter case, the last en- 
tries of the given 4x4 matrix are associated with the p^-like 
orbital. Note that the values correspond to the entries of 
the real-space Wannier (-like) hamiltonian, i.e., a minus sign 
is included. 



chemically meaningful, the orbital filling for this model 
amounts to hve electrons in the present cases, i.e., one 
more than the half-hllcd sp valence of the carbon-group 
elements (Si,Gc,Sn). Since equipped with similar elec- 
tronegativity, a significant charge transfer in the strongly 
covalent surface systems is not expected. Still the five- 
electron filling may be acceptable as the Sn locality of 
the (sp 2 + p z )-like orbitals is neither very strong nor 
demanded on formal grounds. In addition we find this 
model more justified than a simpler three-band model for 
p-likc orbitals only, which would in the end amount also 
to a five-electron filling for three effective bands. Such 
an extreme filling of 5/6 would render the interpretation 
in semi- local terms rather difficult, besides the inconsis- 
tency with the facts of the chemical bonding. 
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B. Three-site Sn unit cell 

It is believed from experimental studies that the struc- 
tural ground state of the Sn/Ge(lll) system corresponds 
to a two-down-one-up (2D-1U) distortion of the Sn sub- 
monolayer—. To be able to account for this reconstruction 
(cf. Fig. [1]), we thus performed electronic-structure cal- 
culations using a 3x3 supercell with a difference A=0.32 
A between the up-down Sn positions^, which we applied 
to our PBE-GGA structurally relaxed unit cell with the 
planar submonolayer. In these extended supercell cal- 
culations only one orbital per Sn adatom is included in 
the subsequent MLWF construction. This results in a 
Kohn-Sham-Wannicr hamiltonian which corresponds to 
a single-orbital problem on a three-site triangular clus- 
ter. For comparison, we derived such enlarged hamilto- 
nians also for the flat systems via corresponding larger 
unit-cell calculations for planar Sn/(Si,Gc)(lll). The re- 
sulting Wannicr-like bands are shown in Fig. [5J Whereas 
the bands again fit exactly in the Sn/Si(lll) case, for the 
distorted Sn/Ge(lll) surface a shift towards lower lying 
bands at the T-point is visible. The latter feature is again 
due to the band hybridizations already observed in the 
single-site Sn unit cell. It may be observed that the 2D- 
1U surface reconstruction leads to a small splitting of the 
low-energy bands, affecting mainly the occupied part of 
these states. 

Besides the study of the influence of apparent sur- 
face reconstructions, the investigation of the magnetic 
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FIG. 8: (Color online) PBE-GGA band structure (dark) with 
extracted Wannier(-like) bands (cyan/grey) for the 3x3 su- 
percell. Top: planar Sn/Si(lll) system. Middle: planar 
Sn/Ge(lll) system. Bottom: Sn/Ge(lll) with the 2D-1U 
reconstruction. 



behavior seems a most important endeavor. As already 
outlined in the introduction, the quasi-2D triangular Sn 
submonolayer might be a realistic case for the appli- 
cation of model ideas discussed in the context of 2D 
quantum magnetism. Therefore we investigated possi- 
ble magnetic orderings within PBE-GGA, especially for 
the Sn/Si(lll) system where intricate ordering patterns 
are heavily discussec&Si. Although (finite-temperature) 
magnetism would be very interesting, one has to keep 
in mind that the constituents (Sn,Si,Ge) arc no high- 
susceptible magnetic materials. The sp-bonding (involv- 
ing large principal quantum numbers) with filled d states 
renders magnetism from a chemical point of view ques- 
tionable. Furthermore the rather large Sn-Sn nearest- 
neighbor distance of ~6.7A asks for a robust exchange 
path to facilitate the appearance of long-range order. 
Still the low-dimensional 2D-charactcr might be sufficient 
to induce delicate magnetic behavior. Note also that the 
high DOS close to the Fermi level (cf. Fig. [5]) might give 
rise to flat-band ferromagnetis m 48 i 49 . 

The results of our investigation of different magnetic 
orderings on planar Sn/Si(lll) 3x3 are summarized in 
Tab. [TTJ The calculations show that especially the fer- 
romagnetic (FM) order is rather intriguing. One may 
stabilize a FM solution for Sn/Si(lll), resembling previ- 
ous work by Profetta and Tosatti^, however this state is 
energetically unfavorable compared to the nonmagnetic 
(NM) solution. The corresponding local energy mini- 
mum of this FM state appears to be rather flat, thus al- 
ready small disturbances drive the DFT self-consistency 
cycle towards the NM state. This was confirmed within 
all three utilized band-structure codes, i.e., MBPP, CP- 
PAW and VASP, whereby the Sn-substrate distance does 
not influence this qualitative result. Note the rather 
small local Mg n moment compared to total FM moment 
of the supercell. Thus this metastable FM state is far 
from being originated from pure local Sn moments, but 
has significant nonlocal character. Since collinear AFM 
order is impossible due to frustration on the undistorted 
lattice, a collinear ferrimagnetic ordering (two up spins 
and one down spin on the minimal triangle) as well as 
the in-plane noncollinear 120°-AFM state were investi- 
gated. Indeed, both latter ordering patterns are found 
to be stable with respect to the nonmagnetic solution, 



magnetic ordering 


E [mcV] 


M [fx B ] 


M Sn [mb] 


ferromagnetic 


1.5 


0.80 


0.031 


collinear ferrimagnetic 


-1.6 


0.44 


0.055 


120° antiferromagnetic 


-3.6 


0.00 


0.058 



TABLE II: Comparison of the different magnetic orderings 
in the flat Sn submonolayer of the V3x\ /r 3R30° surface with 
the Si substrate. The energies E and total moments M are 
given with respect to the nonmagnetic structure in PBE-GGA 
and correspond to an enlarged 3x3 unit cell. This cell incor- 
porates 66 atoms, namely 54 Si, 3 Sn and 9 H atoms (for 
saturation of the bottom bulk-like Si layer). 
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with the lowest total energy for the 120°-AFM ordering. 
Though the local Sn magnetic moments Ms n are only of 
the order of ~0.06 within PBE-GGA, detailed conver- 
gence studies elucidated nonetheless their nonzero value. 
Note that these local Mg n are supplemented by addi- 
tonal spin-polarisation on the remaining sites and the in- 
terstitial contribute. Hence again the picture of strictly 
localized magnetism is not appropriate on the weakly- 
correlated modeling level, but still local moments with 
small AFM exchange may exist in these Sn submonolay- 
crs. Those may possibly give rise to spin-liquid physics or 
eventual magnetic long-range order also in the correlated 
regime^. 

Similar magnetic PBE-GGA studies for a model planar 
Sn/Ge(lll) 3x3 surface did not result in magnetic long- 
range order. The stronger hybridization of the Ge(4s4p) 
states with tin should generally weaken the magnetic ten- 
dencies compared to the silicon substrate. We did not 
investigate the more realistic 2D-1U resonstructed struc- 
ture in this matter, which is left for further studies. 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.: 
CO (eV) 



FIG. 9: (Color online) Realistic single-site DMFT(QMC) 
results for the planar one-band models at /3=200 eV -1 
(T~ 58K). 



IV. INVESTIGATION OF ELECTRONIC 
CORRELATIONS 

In order to take the principal effect of electronic corre- 
lations into account, we concentrate in the following on 
the realistic one-band models derived in section UlI Al and 
their cluster extension from section 1111 Bl Furthermore 
we restrict the investigations to paramagnetic modelings, 
i.e. do not cover possible magnetic orderings. Since it 
became clear from the last section IIII Bl that the energy 
scale for magnetic long-range order is rather small, such 
an analysis shall be postponed to future studies. 



A. LDA+DMFT(QMC) study 

By combining our Wannier hamiltonians with the 
DMFT framework we are in the position to reveal the 
spectral function of the surface systems in the interact- 
ing regime. Due to the small bandwidth of VF~0.4 eV, 
already small absolute values for the Hubbard U may 
introduce strong correlation effects. Although the low- 
energy bands of the discussed systems are mainly com- 
posed of Sn(5s,5p) states, an onsite Coulomb interaction 
of this order of magnitude may very well be reasonable. 
Note that also the screening capabilities because of the 
semiconducting substrate are limited. The results of the 
paramagnetic single-site DMFT(QMC) calculations for 
the realistic one-band model are displayed in Fig. [9] for 
different choices of the Hubbard U. Of course, due to the 
simplicity of the modeling the spectral function (finally 
obtained via the maximum-entropy method) follows the 
usual behavior with increasing U, i.e., a low-energy band 
narrowing with the additional appearance of Hubbard- 
band features at higher energies takes place. The Mott 
transition is reached at a critical value of about £/ r ~0.6 
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Sn/Ge cluster 




0.0eV(d) 


0.819/0.802 


0.0eV(u| 


1.378 


0.5eV(d) 


0.989/0.988 


0.5 eV (u) 


1.022 


0.6eV(d) 


0.996/0.995 


0.6 eV (u) 


1.008 


0.7eV(d) 


0.997/0.997 


0.7 eV (u) 


1.005 



'.8 -0.6 -0.4 
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FIG. 10: (Color online) Realistic cellular-cluster 
DMFT(QMC) results for planar Sn/Si(lll) and 2D-1U 
Sn/Ge(lll). The up/down tags refer to the inequivalent 
atomic positions and the numbers yield the orbital occu- 
pation. Note that albeit slightly different by symmetry, we 
only show the averaged spectral function for the two 'down' 
atoms (occupations are site selective though). Calculations 
were performed at /3=70 eV" 1 (T~ 166K). 



eV for the Sn/Si(lll) surface systems (with a somewhat 
larger value for planar Sn/Ge(lll)). Since this transition 
shows a weak first-order character in the calculations (see 
also section HVBj) critical interactions t7 cl , U C 2 govern 
this regime, respectively. However in the present study, 
due to the simpicity of the one-band modeling, we did 
not map out the respective hysteresis loops. The en- 
ergy gap within the insulating Mott state is of the order 
of ~^0.2 eV for £/~0.7 eV with Hubbard bands at around 
0.3 eV above and below E F for both Sn/(Si,Ge)(lll) sys- 
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terns. The position of our Hubbard excitations is roughly 
in line with the region of increased spectral- weight trans- 
fer measured in photoemission experiments^ below 30K. 
Hence a rather small Hubbard U, slightly larger than the 
single bandwidth, is sufficient to drive the Sn-dominatcd 
surface band Mott insulating. This value is however not 
very suprising since from the involved quantum numbers 
of the respective states no large value of U (i.e. as for 
transition-metal and/or / systems) is expected. 

In order to study the importance of intersite self- 
energy effects, especially for the distorted Sn/Ge(lll) 
case, we also performed cluster DMFT (CDMFT) (for 
recent reviews see e.g. Ref. (50l - l52j ) computations within 
the cellular cluster framework for the basic Sn trian- 
gle in the submonolayer. The resulting spectral func- 
tions arc shown in Fig. [TO] Note that we employed a 
somewhat higher temperature within the QMC solver be- 
cause of the larger numerical effort in the cluster frame- 
work. For the planar Sn/Si(lll) case the main differ- 
ences to the single-site results are given by a slightly 
smaller [/ c ^0.55 cV, smaller energy gap and Hubbard 
bands in some closer range to the low-energy region. We 
checked that those observations are not only due to the 
different temperatures of our computations. One may 
also observe a stronger asymmetry in the spectrum up 
and below the Fermi level, with some stronger reduc- 
tion of low-energy spectral weight below Ep. For the 
site-resolved spectral function of the Sn/Ge(lll) 3x3 
system in the distorted 2D-1U structure it is first im- 
portant to remark that the two downwards shifted Sn 
atoms are inequivalent by symmetry due to different hy- 
bridization with the 'up' atom even assuming the same 
height for the 'down' 2D atom o 19 i 53 (see section HVBI) . 
The latter are moreover less occupied than the upwards 
shifted Sn atom£. With increasing U, this filling inbal- 
ance because of the distortion-induced crystal-field shifts 
becomes smaller and is compensated at the Mott transi- 
tion (see Fig. ITU1) . However a slightly larger critical U c 
compared to the planar cluster for Sn/Ge(lll) is neces- 
sary to reach this transition (more or less equivalent to 
the critical U within the planar single-site DMFT). The 
above noted spectral asymmetry in the occupied and un- 
occupied part is even larger in the distorted case. 



B. LDA+RISB study 

In addition to DMFT(QMC) computations we have 
performed slave-boson calculations within the RISB 
framework to further verify our results and to achieve 
a better resolution of the qualitative differences in val- 
ues for the critical U in the various cases. Furthermore 
we also want to elucidate the intersite spin correlations 
in the cellular cluster modeling. The RISB results for 
the quasiparticle weight Z arc shown in Fig. III! In the 
one-band case, the RISB method yields for both flat sys- 
tems a first-order transition from the paramagnetic metal 
to the paramagnetic insulator. The corresponding crit- 




0.2 0.4 0.6 0.8 

U(eV) 



FIG. 11: (Color online) RISB results for the quasiparticle 
weight Z for the planar single-site (top) and the cellular- 
cluster (bottom) models. 



ical U values are given by J7 C ~0.75 eV (0.78 eV) for 
for Sn/(Si,Ge)(lll) V3x\/3.R30°, hence show the same 
qualitative trend as DMFT(QMC). A somewhat larger 
absolute value within RISB is understandable from the 
saddle-point approximation which is identical to the ne- 
glection of quantum fluctuations. The cellular-cluster in- 
vestigations again verify the reduced critical Hubbard U 
values (similar to the ones from DMFT(QMC)) due to an 
increase in the correlation strength via the inclusion of 
the nearest-neighbor self energies. Also the enhanced U c 
for the distorted 2D-1U structure of Sn/Ge(lll) is a solid 
result. Note that the first-order character of the transi- 
tion is straightened in the cluster modeling. To allow for 
a commensurable Mott transition, the correlations have 
to drive additional charge transfers between the Sn site 
in the 2D-1U structure. The therefore enhanced charge 
fluctuations in the latter case may thus account for the 
larger Uj^. In Fig. [T^] we show the continuous develop- 
ment of the individual site fillings with increasing U for 
the 2D-1U structure. The corresponding labeling of the 
different Sn atoms is given in the inset of Fig. As 
it is clearly seen, the originally enhanced occupation of 
the low-energy orbital for the upwards shifted Sn atom 
on the cluster becomes reduced with increasing U, while 
the downwards shifted ones gain electron filling in their 
respective orbitals. This is in line with the results ob- 
tained from the more elaborate QMC solver to DMFT 
(cf. Fig. [10]). 

In order to provide some insight into the magnetic be- 
havior with taking into account electronic correlations, 
Fig. Q2] displays the local spin correlations (S^Sj) be- 
tween the Sn atoms from the low-energy modeling. It is 
seen that the spin correlations are always negative, i.e., of 
AFM character, as expected by considering the superex- 
changc induced via U. Naturally, the degree of magnetic 
behavior is therewith increased compared to the nonin- 
teracting case. The comparison of the two different sub- 
strates with the planar geometry of the Sn submonolayer 
exhibits the stronger magnetic correlations within the 
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FIG. 12: (Color online) Orbital occupations of the different 
Sn atoms in the 2D-1U structure for Sn/Ge(lll) obtained 
from the cellular cluster calculation with varying U. Inset: 
Designation of the different Sn atoms in the cellular cluster 
description of Sn/Ge(lll) from a top view on the surface. 
The (red/grey) dashed lines mark the cluster utilized for the 
cellular modeling. 

Sn/Si(lll) system. Thus the qualitative result obtained 
from the weakly-correlated PBE-GGA modeling extends 
to the strongly correlated treatment. In the 2D-1U struc- 
ture of Sn/Ge(lll), the values of (S;Sj) for the now 
different Sn-Sn pairs show interesting behavior. Some- 
what counterintuitive, the spin correlations between the 
two up-down pairs, here denoted Sn(l)-Sn(3) and Sn(2)- 
Sn(3), scale rather differently with increasing U. While 
for the Sn(2)-Sn(3) pair (S^Sj) is most negative, for the 
Sn(l)-Sn(3) pair the spin correlations become nonmono- 
tonic close to U c , anticipating tendencies to eventual FM 
coupling. The remaining Sn(l)-Sn(2) pair of the two 
down atoms scales inbetween these functions. This result 
can be explained by the fact that the absolute value of 
the hopping amplitude t is maximum between the Sn(2)- 
Sn(3) pair and is accordingly weakened between Sn(l)- 
Sn(3). Because of the upward shift of the Sn(3) atom, 
the lobe from its effective Wannicr orbital pointing to- 
wards Sn(2) achieves a wider range, leading to stronger 
overlap with the Wannier orbital of Sn(2). 

V. SUMMARY 

We have presented a detailed modeling of the 
Sn/(Si,Ge)(lll) surface systems using state-of-the-art 
band-structure methods combined with many-body tech- 
niques. In addition to previous DFT investigations 
we extracted realistic one- and four-band Wannier (-like) 
hamiltonians that are capable of describing the electronic 
structure close to the Fermi level. The physics of the 
minimal one-band model is of course restricted, how- 
ever due to the prominent Sn-dominatcd half-filled sur- 
face band at Ef, especially for Sn/Si(lll) y/3x^/3R30° 
this is believed to be an adequate modeling viewpoint. 
The atomistic description of the Sn/Ge(lll) system still 
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FIG. 13: (Color online) Spin correlation functions for the 
Sn/Si and the Sn/Ge distorted cluster obtained from the 
RISB study. Top: planar systems, bottom: distorted 2D- 
1U Sn/Ge(lll) system. The vertical line marks the critical 
U value, respectively. 

raises some questions concerning the apparent limitations 
of simplified exchange-correlation functionals in this case. 
Nonetheless on the present level of investigation and com- 
parison, the PBE-GGA perspective appears sufficient to 
reveal the essential differences between both surface sys- 
tems. The hybridization of the Sn submonolayer with 
the substrate is stronger in the case of Ge, leading there 
to a more intriguing entanglement of that layer with the 
supporting atoms below. Consequently, structural re- 
constructions are more likely and indeed take place as 
verified in experimental studies. 

The magnetic behavior in these systems renders them 
most fascinating from a fundamental physics point of 
view in terms of possible low-dimensional quantum mag- 
netism. The (non)collincar spin-polarized DFT calcula- 
tions revealed that there is indeed the chance for intricate 
magnetic ordcrings. In this respect we found the 120°- 
AFM structure to be the magnetic ground state in the 
weakly-correlated description. Further theoretical stud- 
ies, including the effect of true many-body correlations on 
the magnetism are important and should be motivated by 
the present work. Although the resulting moments are 
rather small, maybe advanced experimental techniques 
are capable of exploring these small energy scales. 

Using DMFT(QMC) and RISB techniques, a moder- 
ate value of the order of U ^0.5-0.6 eV was found to gov- 
ern the systems from a one-band Hubbard-like model in 
a strongly-correlated limit. Because of the small low- 
energy bandwidth it is reasonable to expect electronic 
correlations to be important, however sole local Coulomb 
interactions are surely limited in the description of these 
systems based on carbon-group elements. Yet already 
on this rather simple level a delicate interplay between 
the structural and electronic degrees of freedom, antici- 
pated in other model(-like) framcwork o 13 ' 15 , was verified 
on the grounds of an LDA+DMFT approach. Future 
studies now have to concentrate on including the sub- 
strate degrees of freedom explicitly in the modeling of 
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the interacting electronic structure, perhaps even allow- 
ing for the dynamical structural reconstructions needed 
for Sn/Ge(lll). 
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